NASA Contractor Report 141439 


The Prediction and Mapping of Geoidai 
Undulations from GEOS-3 Altimetry 


William Kearsley 


(lJaSA-CE-141439) THE PREDICTION AND HAPPING 
OF GEOIDAI DNDOIATIONS FROM GEOS-3. AITIHETRI 
(Ohio State Oniv. Research Foundation) 52 p 
HC A04/HF A01 CSCI 08B 


N78-2'-jJ07 

i 


Dnclas 

G3/43 ; 16685 


April 1978 


IMASA 

Nationai Aeronautics and 
Space Adnninistration 

Wallops Flight Center 

Wallops Island, Virginia 23337 
AC 804 824-3411 




NASA Contractor Report 141439 


The Prediction and Mapping 
Undulations from GEOS-3 Altimetry 


William Kearsley 
The Ohio State University 
Research Foundation 
Columbus, Ohio 4321 2 


Prepared Under Contract No. NAS6-2484 


rvv\s/\ 

National Aeronautics and 
Space Administration 

Wallops Flight Center 

Wallops Island, Virginia 23337 
AC 804 824-3411 


Geoidal 



Foreword 


This report was prepared by Dr. William Kearsley, Research 
Associate, Department of Geodetic Science, The Ohio State University, under 
NASA Contract No. NAS 6-2484, The Ohio State University Research Foundation 
Project No. 783904, Project Supervisor, Richard H. Rapp. The contract 
covering this research is administered throu^ the NASA Wallops Flight Center, 
Wallops Island, Virginia 23337, with Mr, Ray Stanley, Project Scientist. . 


-ii- 



Table of Contents 


Foreword , ii 

Acknowledgements iv 

1. Introduction 1 

1.1 Map Specifications 1 

1.2 Description of Programs 1 

2. Filtering the Altimetiy Data 2 

3. Interpolation of Data for Contouring 4 

3.1 Least Squares Prediction 4 

3.2 Prediction by Weighted Means 5 

4. Choosing the Grid fiiterval for Contourii^ 16 

5. The Predicted Geoid - Results and Accuracies 18 

5. 1 Accuracy of the Geoid from Altimetry 18 

5. 2 Results of Computations 19 

6. Conclusions 62 

Bibliography 63 

Appendix A - Instructions for Use of Contouring Program 65 


Appendix B - Flow Chart - Main Program 


68 



Acknowledgements 


I am indebted to Professor Eichard H. Rapp for the opportunity to 
work with him in the analysis of the GEOS-3 data, and also for his help and 
encouragement on this project. I would also like to thank Laura Reiter, for her 
care and patience in typing fiie manuscript of this report. 


-iv- 



1. Introduction 


This paper is written to describe tiie methods developed at the 
Department of Geodetic Science to produce 'geo id' maps over ocean areas using 
the adjusted altimetry data described in Happ (1977). This data represents the 
results of bias and trend removal from the original, but edited, Geos-3 altimeter 
data supplied by NASA through the Wallops Flight Center. 


1. 1 Map Specifications 

The data forming the basis of the contour maps is the estimated 
geometric sea surface height (h) which contains corrections for tropospheric refraction 
and altitude bias. This provides an estimate of the separation between the sea 
surface and the actopted reference ellipsoid. IE the effects of ocean tides can now 
be removed from h, and medium and short wavelengQi noise filtered, the result 
will be an estimate of the geoid-spheroid separation, N. To find the best contour 
interval to adopt for a contour map of the geo id, one must first estimate the error 
likely to be present in N (^N). 

The standard deviations of the adjusted altimeter data are due to 
instrumental effects and the unmodeled orthometric he^ht, and are estimated to be 
about ±1 m for short pulse and ±2 m for long pulse measurements (Rummel, 1976, 
p. 12-13). These accuracies can be used to determine a suitable contour interval 
for geoidal maps based on the altimeter data. For example, the contour interval 
is usually 3. 33 6 n (American standards to 5 6 n (European standards), 6 n being 
the error in individual heights (Richardus, 1973, p. 83), suggesting a contour 
interval of 3 m to 5 m for data measured in the short pulse mode. However, these 
standards are for medium to large scale maps. Generally, the geoid is smooth 
and geoid maps are small scale (e. g. 1:500,000), so it appears reasonable to adopt 
a 2 m contour interval in this case. This compares well with (Ibid. , p. 106) 
where 2 6 n appears to be sufficient contour interval in flat areas. 


1. 2 Description of Programs 

The programs used in the prediction and contouring of the undulations 
are referred to throughout the report and are listed here with a brief description 
for convenience. They were written for the IBM 370 at the Ohio State University, 
which is ecpipped with an on-line VERSATEC plotter. The details for user operation 
are given in Appendix A, and a Flow Chart of the program MAIN can be found in 
Appendix B. 
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Programs used in the production version of the contouring package include; 


(i) MAIN which reads all altimetry data from the relevant tape, and having passed 

.this data through a window specifyir^ the map limits, organizes the 
subroutine used to predict geo id values, plot the contours and axes, and 
so on. (See Appendix B for Flow Chart) 

(ii) GPP the general plotting package developed by the Computer Services Unit to 

facilitate plotting on the VERSATEC plotter. 

(iii) PRDT a subroutine originally developed by Snowden (1969), but largely 

rewritten to suit the present task and system more closely. 

(iv) GEOGRD- a subroutine to convert geographical to projection coordinates using a . 

spherical model of the earth of radius 6.371 x 10® m. Options include the 
Plate-Carree, Mercator, Cassini-Soldner and Equal Area Conical and 
Cylindrical projections (see also Pearson, 1977). 

(V) COVA a subroutine to compute global or local covariance functions (Tscheming 
and Rapp, 1973). 

(Vi) CGTOUR a subroutine developed by Snowden (1969) which plots and annotates 
contours from gridded data. 

(vii) AXPLT a subroutine developed to plot and annotate rectangular map borders. 

(viii) SHIN a subroutine developed to plot the coastlines of those land masses 
falling within the mapping limits .. 


2. Filtering the Altimetry Data 

The data to be used in this report has been processed to minimize the 
large long wavelength and constant biases in the original data. The information which 
is now relevant to Ihe contouring process is (i) the measurement type, (ii) the sequence 
number, (iii) geocentric latitude, (iv) geocentric longitude, (v) the estimated standard 
deviation of the range and (vi) the adjusted sea surface height (for these purposes, 
equivalent to the geoid height, N), 

The adjusted data contains high-frequency noise which should be filtered 
before it is used in production (Rummel, 1976, pp. 12-15). A least-squares filtering 
technique developed by Rummel was used for this purpose. Using a band width of 
11 points the filter is passed through each segment of an orbit, such that data from 
any adjacent segment is not used in the filtering process." 
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Table 1 shows typical filters developed for an 11 point band width. 

Filter No, 1 is applied when the subject data is the first in the 11 point section. 

Filter No. 6 where the subject data is the sixth, and so on in reverse order (e. g. 
Filter 5 for point 7). The filter values are generated from subroutine COVA 
(Section 1. 3) based on the lag of the data points averaged over the segment for 
test purposes. 

The data for the Alaska test area, which resulted from a preliminary 
adjustment, was also passed through a filter usit^ a moving averages technique 
and the results compared with those from the least-squares method. Using a 
sample of about 400 points it showed that 85% of the values from the moving averages 
method differed by less than 0,1m, and 95% differed by less than 0,15 m. The 
largest differences (0,25) occurred at the terminals of the sections. This suggests 
that the moving-averages method introduces no significant error (cf. the least 
squares approach) in data being used to plot 2 m contours (see Section 1.2), and 
costs about 2/3 the cost of the former approach. 


Table 1 

Coefficients of Typical Filters Generated from COVA 
for Least Squares Filtering 


Filter 

# 

Filter Values 

r— 1 

0.24 

0,20 

0.16 

0.13 

0.10 

0.07 

0.05 

0.03 

0.01 

0.00 

-0.02 

2 

0.20 

0,18 

0.16 

0.13 

0.10 

0.08 

0.06 

0.04 

0.02 

0.01 

0.00 

3 

0.16 

0.16 

0.15 

0.13 

0.11 

0.09 

0.07 

0.05 

0.04 

0.02 

0.01 

4 

0.13 

0.13 

0.13 

0.13 

0.11 

0.10 

0.08 

0.06 

0.05 

0.04 

0.03 

5 

0.10 

0.10 

0.11 

0.11 

0.12 

0.11 

0.09 

0.08 

0.07 

0.06 

0.05 

6 

0.07 

0.08 

0.09 

0.10 

0.11 

0.11 

0.11 

0.10 

0.09 

0.08 

0.07 
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A question arises which concerns the order in which the adjustment and 
the filterir^ is performed. The long wavelength errors (n<20, where n denotes 
the degree of a Legendre polynomial ejqjansion between - tt and + tt) result from 
errors in the computed orbits. The medium (20 ^n ^ 200) and short (n > 200) wave- 
length errors are due mainly to sea surfe.ee state and altimeter noise. In theory 
one should remove these errors before the external constrains at crossover points 
are placed on the data during the adjustment. In practice, it was found too expensive 
to filter the data, and the need for such a filter was diminished by the fact that 
(i) the data had already been passed through a coarse filter (Rapp, 1977, pp. 1-, 2) 
and (ii) the process of predicting the grid point value itseK performs a smoothing 
or filtering function. 

As a result, all maps produced once tests had been completed and the 
program finalized (i.e. Figures 6 - 45) used adjusted data which had not been 
fdtered. 


3. Interpolation of Data for Contouring 
3. 1 Least Squares Prediction 


The method of observation results in the data being densely spaced along 
the satellite track (0. 3 degree), with the spacing between tracks in the test area 
being comparatively large and irregular. Attempts to contour direct from the data 
feiled, and it was necessary to interpolate the data onto a grid pattern in order to 
t£dce advantage of the subroutine CGTOUR developed by Snowden (1969). 

The subroutine PRDT was developed for this purpose. This generates 
a grid, the dimensions of which are specified by the user and then, using a least 
squares prediction technique, estimates a value for the geoid height at each of the 
grid points. An estimate of accuracy is also given for this value. This information 
can now be used by subroutine CGTOUR to obtain and plot geoidal contours for the 
area. 


A note should be made here concerning the accuracy estimates. These 
will be based on (i) the standard deviation of the data used in the prediction and (ii) 
the separation of these data points from the prediction (grid) point. The standard 
deviations in (i) are supplied with the original data and are a measure of internal 
consistency (rather than of absolute accuracy) of the repeated observations used to 
obtain the mean altimeter range value. The prediction error in (ii) is a test of how 
well the data fits the statistical model assumed in the least squares prediction. The 
accuracies computed can not therefore claim to be absolute but are a good indication 
of the precision or relative accuracies of the predicted values. The value of having 
this information can be appreciated by referring to Figure 2, where one can clearly 
see how the errors are related to the tracks of the satellite orbits (and hence the 


- 4 - 



separation between the prediction and data points) . For comments relating to the 
accuracy of the contour map, see Section 4, The results of the predictions and 
their errors along with their respective plots, for a 15° x 10° region in the Alaska 
test area, are shown in Tables 2 and 3 and Figures 1 and 2. 

The basis of the least squares prediction is the covariance fimction 
which expresses the covariance relationships between any two geoid undulation 
values. The subroutine COVA (Tscheming and Rapp, 1975) was used to generate 
these covariance values . Attempts were made to simplify and speed up the 
computation by using simple polynomial models such as the second and third order 
Markov models suggested by Jordan (1972), and the Gaussian model by Moritz (1972), 
but the time saving was negligible and in some cases, the predicted values were 
obviously wrong. For example, in one instance the value differed by about 300 m 
from that predicted by COVA. This erratic behavior was caused by instability of 
the auto-covariance matrix generated by the third-order Markov function, which 
is uncritical for sn^ll changes in 0 as 0-*O. The second-order Markov function 
is more sensitive in this region and was found preferable in practical computation, 
even though it is felt to be a poorer reflection of the true behavior of the covariance 
function for small 0 (Jordan, 1972). 

Another alternative has been suggested by Lachapelle (1976, p. 7 ), who 
found a saving of a factor of 2. 5 to 3 when a tabulation of covariance values generated 
from COVA was used. When this was adopted, some saving was noted, although not 
nearly as dramatic as that mentioned above. The apparent reason for this is that 
only the total cost of the program is known, and part of it-that of contouring the 
gridded geoid height and accuracy values-remains constant and fairly expensive. 

The final version of the program uses a tabulation of COVA values at every 0.05 
degrees, 0 ^0 ^ 10, and employs a linear interpolation to extract C (0) when 0 is 
not an even 0?05 value. 


3. 2 "Prediction by Weighted Means 

The simple deterministic model: 



1=1 
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where 


A 

Nj = 
Ni = 
Tij = 

u = 
and n = 


predicted value of N at point j 
value of N at point i 
separation between points i and j 
power of prediction 
number of observations 


has been used by Sjoberg (1975, p. 7) to predict gravity anomalies (also see Davis, 
1973, pp. 316-318) and it was decided to test it in the estimation of geo idal undula- 
tions at the grid points. Sjoberg shows that the optimum value for u is about 3. 5 
(Ibid.), and tests were performed with the Alaska Data Set to see what effect the use 
of this simple model would have on the resultant map. Two values of u were tested 
(i^= 1 and l^= 3. 5) for both a 5 data points and a 10 data points prediction. It was 
found that, in this case, the value of v had little influence on the resulting prediction. 
The use of 5 data points in the prediction was about 20% less expensive than 10 data 
points (whose cost was comparable to that of the least squares prediction). However, 
there appeared to be little advantage in using 10 points in the prediction as the 
resulting values were significantly unchanged. 

An example of the maps produced by the simple wei^ted mean method 
is given in Figure 3. 


Errors in the Predicted Quantities 


The error in the predicted value (CTr) will come from two sources, (i) the 
error in tiie data (Od) and (ii) the error of representation of the weighted mean (Cp). 
These two errors are uncorrelated for all practical purposes, and the total prediction 
error is expressed as: 


(3.1) 


3 

Or = 


+ Ud 


CTd 

value from the 


can be gauged from an analysis of the differences of each predicted 
mean value . Thus : 


(3. 2) 


3 

(Td = 




pvv 
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where 


i = 
and p 1 = 

Ti = 


1 to no. of points used in prediction 

1/r,^ 

separation of point i from the prediction point. 


CTp is related to either the size of the area for which ihe predicted 
value represents a mean value, or the 'distances' over which the prediction is 
made. In both cases, a function of distance (d') is involved, and thus we can 
say: 


3 

CTp = 


c d' 


For convenience in this case (since this quantity was already computed) 

(3.3) d' = ^ rf" 

1 

and by experiment c was found to be. about 30 m®. 

The final result 

^ p V V 

= i + 30 

I p(i-i) 

was used to estimate the error of prediction at each grid point and the resultant map 
of the errors showed the same, trends as for the least squares error analyses (compare 
Tables 3 and 5, and Figures 2 and 4). Table 5 and Figure 4 are derived from 10 data 
points in the prediction of each grid value. The errors estimated using 5 data points 
are considerably higher, as might be expected, although the actual estimates of N 
are significantly the same. 

The differences between the least squares prediction and the weighted 
mean (with 10 points) gives an RMS of about 1.5 m. In areas where the data is 
comparatively dense (see Figure 1) the differences are small (compare,' for example, 
rows 7 and. 8 in Table 4 and Table 2), but in some instances they are unacceptably 
large. For example, point (16, 8) is 0.0 in Table 2 and -13.8 in Table 4. This 
occurs, in an area of no data and reflects the weakness of the methods when applied 
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Figure 1 
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Figure 3 


Geotdal Undulations For Alaska Test Area Using Weighted Means 

Contour Interval: 5 m 
Projection; Cassinl-Soldner 
Grid; 21 x 11 
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Table of Predicted Geoid Heights by Least Squares Prediction 
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Table 3 


Table of Predicted Standard Errors 
Alaska Test Area 
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Table 4 


Table of Predicted Geoid Heights by Weighted Means: 

u = I 

(lO Closest Points) 
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Table 5 


Table of Predicted Standard Errors, by Weighted Means 

v= 1 ■ 

(Using 10 Closest Points) 
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to a situation requiring extrapolation. It appears that, providing the area is densely 
covered by data, a simple we^hted mean approach produces quite adequate results. 
However, the total cost of the prediction are similar, and the advantage of using 
the least squares method is that it provides a super-structure of theory which is 
lacking in the more intuitive weighted means approach and also, that it simplifies 
the computation of liie expected error of the prediction. 


4. Choosing the Grid Interval for Contouring 

The usual specifications for a medium to large contour map is ihat contour 
interval (I) bebetween3»33 Oh and 5 Oh , where CTh is the standard deviation of the 
individual heights in the mappii^ area (Richardus, 1973, p. 83). While the maps 
involved here can hardly be considered medium or large scale, it is considered 
reasonable to adopt as the standard I = 20 h > as the comparative smoothness of the 
geo id over large areas makes it unfeasible to apply small scale standards to this 
problem. The internal consistency of the altimetry data gauged from the discrep- 
ancies in the crossover points after adjustment (see Rapp, 1977) suggests a Oh of 
about 0. 5 m,but it appears that a contour interval of less than 2 m would be unrealistic. 
In fact, in earlier tests a contour interval of 5 m was used, but with the refinement 
of the technique, 2 m was finally adopted. 

In normal topographic mapping the contours are either interpolated 
linearly between points defined by field survey at changes of grade, or are plotted 
directly from a three dimensional photographic model. The criteria applied in these 
two cases are not helpfol in this present instance, where the contouring will be done 
from height values predicted on a pre-determined grid pattern. The problem is to 
find the optimum grid pattern which will; 

(i) . faithfully represent the mapping area in a discrete form 

(ii) produce realistic contour lines 

(iii) not be too costly. 

The first point is probably the most important. Too large a grid interval 
will leave large areas of the map unrepresented. The subsequent contour lines, being 
a simple linear interpolation between grid points will therefore miss the features 
present in the ungridded data, since the grid points lack the important property of 
being located at salient points which adequately describe the features of the geoid. 

The extreme is to make the grid spacing equal to the separation of the altimeter data 
along any one track (=-0. 3 degrees) but such an intensive grid would be very costly 
to predict and map and is unnecessarily dei^e. 
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The normhl criterion for contour line accuracy is that it be accurate 
to within, ± 1/3 i, I being the contour interval. An accuracy of 0. 7 m is bordering 
on the noise level of the altimeter data, and it seems more realistic to prescribe 
an accuracy of 1 m for the contour line. This means that a grid interval should be 
chosen so that it is less than half the wavelength of geoidal undulations having 
amplitudes of at least 1 m (Nettleton, 1976, p. 159). 

This wavelength will obviously vary from area to area, .and it is 
difficult to determine its magnitude without some a priori knowledge of the areas. 

One could use existing geoid maps for ocean regions such as the U. S. Calibration 
Area, but these maps themselves will be limited in accuracy by the grid interval 
chosen as the basis for the map. A study of the adjusted altimetry may reveal 
some of the waveler^h characteristics of the geoid in the mapping area. Failing 
all else, the grid interval can be determined by trial and error. 

Tests with various grid intervals showed a value of 2° failed to sense 
high frequency features related to the ridge, and trench systems, on the ocean floor. 

A 1® grid appeared to be adequate to represent these features and did not increase 
the cost of the maps beyond reasonable limits, m any case, a l° spacing seems 
to be the smallest that could reasonably be used, as spectral analysis of altimeter 
data suggests that below a 100 km half wavelength, the sea surface topography 
seems certain to be dominated by sea-state departures (Wagner, 1977, p. 35). 

With all cost reducing options operating, the cost of producing a 30° x 30° geoid 
map and its associated error map from about 3500 altimeter data was =* $15. 

Some comments concerning the maps of precision are warranted here. 

It must be recognized that these are derived from the expected error of prediction 
(Moritz, 1972, Eq. 4-28a) at the grid points. They are thus maps of the precision 
of prediction at the grid points and may not reflect the actual precision of values 
interpolated between grid points (particularly if the grid interval is too large to sense 
the high frequency features) . Also, they do not reflect the accuracy of the contours 
of the geoid map. This, if the grid interval has been properly chosen and the data 
coverage is adequate, should be of the order of ±lm, as explained above. 

The accuracy of the contours will, in the final analysis, be dependent 
upon the algorithms used in the computer program to interpolate the contour 
lines between the grid points. The routine used in this case (Subroutine CGTOUR, 
Snowden, 1969) joins interpolated points by a series of short straight lines. This 
gave rise to a peculiar 'stepping' or 'staircasing' effect in areas where the grid 
points lay between two widely (i. e. 2°) separated arcs. This was because grid 
points close to one arc used only points from that arc to determine its value, 
while those closer to the second used only data from that arc. In an attempt to . 
overcome this problem, a condition was imposed whereby of the, five close points 
chosen to provide data for prediction, only the first three could come from the ‘ 
closest arc and the remainder had to be chosen from other arcs. 
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5. The Predicted Geoid “ Results and Accuracies 


The errors computed in Sections 3. 1 and 3.2 estimate the precision 
or the internal accuracy of the prediction. To obtain some estimate of the. absolute 
accuracies of prediction, a comparison of the predicted geoid undulations at 
intervals was made with existing geoid values in (i) the U. S. Calibration Area, 
and (ii) the Philippines region of the Pacific Ocean. The results of the compar- 
ison are summarized below. 


5. 1 Accuracy of the Geoid from Altimetry 
(i) U. S. Calibration Area 

The altimeter geoid has been compared with two different geoids in the 
area bounded by 40® 20° and 297° ^ 277°. The geoids used for comparison 

were (i) that described by Rapp and Rummel (1975), where the GEM 6 coefficients 
were combined with the 1° x 1° anomaly data and (ii) the geoid computed by Marsh 
and Chang (1976) who combined the potential coefficients of GEM 8 with l°x 1°, 

15' X 15' and 5' X 5' mean anomalies. 


Table 6 

Undulation Comparisons in GEOS-3 Calibration Area 



Nalt - N?,r 

—1 

1 

O 

Mean Difference 

-0.3 m 

-1.0 m 

RMS Difference 

±1.2m 

±1, 5 m 

[Variancejt 

±1.2m 

±1.1 m 

• Number 

337 

335 

Maximum Difference 

5.0 m 

9.2m 

Minimum Difference 

-3,5 m 

-4.6 m 


The values in Table 6 should be viewed in ihe light of the expected 
accuracies of the Nr.r and Nh,c values. The standard deviations of Nr.r range 
from 1. Im to 2.1.m, with a mean of about 1.6 m, so much of the contribution 
in the differences will come from the Nr^r value. A similar comment can be 
made concerning the Nh, c values. 
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A plot of the differences shows a systematic trend or slope exists 
between the altimetry geoid and the Rapp, Rummel and Marsh, Chang geoids. 

This could indicate small errors in the potential coefficient sets used in the 
analysis (Rapp, 1977, p. 13). However, the comparisons seem remarkable 
considering. the fact that the Nalt contains instrumental noise, and is in reality 
a map of sea surface topography rather than the geoid. The comparison indicates 
that the original estimates of accuracy and the plotting parameters based thereon 
(Sections 1. 1 and 4) are acceptable for continued mapping. 


(ii) Philippines Area (see Figure 6) 

Differences were taken between the altimetry undulation values and those 
supplied by Marsh (1977, private communication) at l°x 1° intersections in the 
Philippines area. Analysis of 625 differences showed ihat the average value was 
-1. 33 m (Nalt - Nharsh)» with an RMS value of ±4.0 m. Upon removing the mean 
difference from these differences, a standard deviation of ± 3. 75 m was obtained. 
Maximum positive difference was 7.4 m while maximum negative difference 
was -13. 3 m. 

The area has a number of high-amplitude features of relatively short 
wavelength (=-3°) and it is unlikely that these features would, be reflected in the 
gravimetrically derived geoid of Marsh. This helps to ej 5 >lain why the comparison 
• in this region is not as favorable as that in the U. S. Calibration Area. 


5. 2 Results of Computations 

The techniques developed above were now applied to adjusted altimetry 
data in the Pacific, Indian and Atlantic Oceans. A basic area of 30°x 30° was 
chosen for mapping (see Figure 5 and Table 7) and the geoid maps themselves, 
with their associated precision maps, are presented in Figures 6' to 45. It is 
felt that for most purposes these maps should provide good service for years to come. 
Areas of the oceans not portrayed in these individual maps did not have sufficient 
altimetry data available at the time of the report preparation to warrant the 
production of the contour maps. 
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Index of Geoid Maps from Geos-3 Altimetry 
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Table 7 


Index - Geo id Maps from Altimetry 


Map 

No. 

Fig. 

No. 

Geog. Limits 

Description 

Lat. 

(degrees) 

Long. 

(degrees) 

I 

6 

30 - 0 

150-120 

Philippines - W. Pacific 

II 

8 


150-120 

New Guinea - N. Australia 


10 

KEEEfl 

297-277 

U. S. Calibration Area 


12 

0 -*-30, 

120- 90 

Indonesia-East Indian Ocean 

EHH 

14 


90- 60 

Central Indian Ocean 


16 


90- 60 

India - N. Central Indian Ocean 

ncH 

18 

30 -* 15, 

280-260 

S. W. U. S. Calibration Area 

musm 

20 

-30^-60, 

120- 90 

S. E. Indian Ocean - S. W. Australia 

IX 

22 

-30 -*^60, 

90- 60 

S. Central Indian Ocean 

X 

24 

-30 -*-60, 

60- 30 

S. W. Indian Ocean 

XI 


O' - -30, 

60- 30 

East Africa - W. Indian Ocean 

xn 

28 

-30 -*-60, 

150-120 

South Australia Basin 

xni 

30 

^gTqBH3ill 

180-150 

New Zealand - Tasman Sea 

XIV 

32 

0 --30, 

180-148 

N. E. Australia - Coral Sea 

XV 

34 

0 -*-30, 

210-180 

Mid Central pacific Ocean 

XVI 

36 

-30 -*-60, 

210-180 

S. Central pacific Ocean 

XXI 

38 

-30 --60, 

330-300 

S. West Atlantic 

XXIV 

40 

MKiggSII 

300-270 

South America 

XXV 

42 

60 -* 30, 

240-210 

N. E. Pacific 


44 

40 - 10, 

260-230 

N. E. Pacific 
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Figure 9 


Geold Undulations in the Philippines 
From Adjusted Altimeter Data 
Contour Interval; 2' meters 
Grid hitarval 1” 

Scale 1; 500,000 
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Figure 7 


Geoid Undulation Accuracy In the Philippines 
From Adjusted Altimeter Data- 
Contour Interval:. 5 meters 
Grid Interval 1* 

Scale 1: 500,000 


LONGITUDE (DEG) 



23 


ORIGINAE PAGE ISf 
OP POOR QUAIOT 




Figure 8 


Geold Undulations - New Guinea 

Grid Interval 1° 

Contour Siterval 2 m 
Scale 1; 500,000 
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Figure 9 


Precision of Prediction - New Guinea 

Grid Interval 1° 

Contour Interval 0. 5 m 
Scale 1; 5,000 
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FlgttW 10. Geoid Undulations in Geos-3 Calibration Area 
From Adjusted Altimeter Data. 

(Contour Interval; 2 meters) 

Grid loteml 1° 

Seals 1: 333,333 
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Figure II. Geoid Undulatioa Accuracy in Geos-3 Calibration Area 
From Adjusted Altimeter Data. 

(Contour Interval: 0.5 meters) 
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Figure 12 

Geoid Undulations - Fast Lidlan Ocean 


Grid Interval 1® 
Contour Interval 2 m 
Scale 1: 500,000 
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Figure 13 

Pieclsioa of Predlcttou ^ East Pu ll a n Ocean 


Grid Interval 1® 
Contour Intervid 0.5 m 
Scale 1: 500,000 
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Figure 15 


Precision of Prediction - Central Indian Ocean 

Grid Interval 1“ 

Contour Interval 0.5 m 
Scale 1; 500,000 
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Figure 16. Geoid Undulations in the Central Sidian Ocean 
From Adjusted Altimeter Data. 

Contour Interval: 2 meters 
Grid Interval 
Scale 1; 500.000 
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Figure 17. Geoid Undulation Accuracy In the Central Indian Ocean 
Prom Adjusted Altimeter Data. 

Contour Interval; 0.5 meter 
Grid Interval l® 

Scale 1: 500,000 
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Figure 18 


Geoid Undulations - S. W. U. S, Calibration Area 

Grid Interval 1° 

Contour Interval 2 m 
Scale 1; 500,000 
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Figure 19 


Precision of Prediction - S. W. U. S. Calibration Area 

Grid Interval 1° 

Contour Interval 0, 5 m 
Scale 1: 500,000 
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Figure 20 


Geold Undulations - S. W. Australia 

Grid Interval 1° 

Contour Interval 2 m' 

Scale 1:500, 000 
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Figure 22 


Geoldal Undulatioiis - S. Central ^dian Ocean 
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Figure 23 


Precis ion of Prediction - S. Central TntHan Ocean 
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Figure 24 


Geoldal TJudulatioiis - S. W. Indian Ocean 

GridSitorval 
Contour Uterval 3-m 
Scale 1: 500,000 
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Figure 25 

PsreciBion of Prediction - S. W. Indian Ocean 


Grid atorral 1° 
Contour Interval 0.5 m 
Scale 1: 500,000 
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Figure 26 


Geoldal UndulatioDS - Zast Africa 
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Figure 27 

Freoislon of Prediction - Fast Africa 
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Geoldal Undulations - Southern Australia 

Grid Qiteival 1° 

Contour Interval 2 m . 

Scale 1: 500,000 
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Figure 30 

Geoldal UnduUtlonB - Tasman Sea 

Grid Interval 1“ 

.. Contour Interval 2 m 
So^tt Is 500(000 
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Figure 32 


Geoidal Undulations - N. E. Australia 

Grid Interval 1“ 

Contour Interv^ 2 m 
Scale 1: 500,000 
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Figure 33 


Precisioa of Prediction - N. E. Australia 


Grid Interval 1° 
Contour Siterval 0.5 m 
Scale 1: 500,000 
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Geoidal Undulations Mid Central PaclBo 
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Figure 35 


Precision of Prediction - Mid Central Pacific 

Grid Interval 1*’ 

Contour fiiterval 0.5 m 
Scale 1: 300,000 
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Figure 36 


Geoldal UndulatlonB - S. Cential Pacific 


Grid mterval 1° 

■ Contour Interval 2 m 
Scale 1: 500,000 
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Figure 37 
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Precision of Prediction - S. Central Pacific 
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Contour Interval 0.5 m 
Scale 1; 500,000 
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Geoidal Undulations - S. W. Atlantic 
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Scale 1; 500,000 
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Figure 41 
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Figure 42 


Geoidal Undulations - N. E. Pacific 

Grid Interval l°x 1° 

■ Contour Interval 2 m 
Scale 1; 500,000 
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Figure 43 

Frecisloa of Prediction - N. F. Pacific 
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Figure 44 


Geoidal Undulatioos - N. E. Pacific 


Grid Interval l“x 1° 
Contour interval 2 m 
Scale 1: 500,000 
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Precision of Prediction - N. B. Pacific 


Grid Ihtarval l“x 1°. 
Contour Merval 0.5 m 
Scale 1; 500,000 
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6, Conclusions 


It is apparent from the foregoing that it is possible to obtain values of 
the geoidal undulation {or its approximation) from GEOS -3 altimetry with an accuracy 
which is of the order of ± 0, 5 m to ± 1. 0 m. It is difficult to ascertain an exact 
accuracy for this parameter as the geo ids used for comparison themselves have 
uncertainties which are- of this same order-of magnitude. 

It also appears that contouring the geoid over the ocean based on 
undulations predicted at 1° intervals produces maps which are representative of 
the geoid. The contour lines are expected to have an accuracy of ± 1 m, assuming 
the grid values themselves are error free. The precision of prediction for each 
geoid map is shown graphically in an associated map, and it would be possible to 
obtain some estimate of the absolute accuracy of the contour lines by combining 
the variance of the prediction with the variance of contouring. In areas well 
covered with altimetry data, and this is the case in most areas included in this 
report, the contouring error is theoretically estimated to be ± 1. 2 m. 

The altimetry data contains signals of higher frequency than the 1° 
half wave-length frequency inferred by the grid spacing. Areas with high frequency 
features such as the Phil^pines Trench or the Tonga Trench would benefit by use 
of a smaller grid interval, perhaps as small as the 0? 3 which approximates the 
spacing of the altimeter data. 

Discrepancies between orbits at crossover points have not (apparently) 
produced any problems in the contourii^ procedure. It is probable that the 2 m 
contour interval is too coarse to show any. effect which.may be produced by these 
discrepancies. 

hi summary, the contouring program produces reliable regional geoid 
maps from the GEOS-3 altimetry data. It is expected that, as more data becomes 
available in areas at present poorly covered, the mapping program will continue 
until all ocean areas of the world will be covered by the standard 30° x 30° geoid 
maps. 
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Appendix A 


Use of the Contouring Program 


A. 1. Data Input 

The altimetry data is read directly Trom the relevant tape, which, in 
the case of the areas mapped in this report, is either GS 102 or GS 110. 

Card 1. Format: 5F10.2, 315, F5.1, 215 

(i) PHIl, PHI2, LAMI, LAM2. (REAL*4) Units - decimal degrees. 

Sets the upper and lower limits in latitude and longitude for the map. 
Also defines the limits of the inner window. Maximum number of points currently 
acceptable is 7,000, 

(ii) OVLAP (REAL*4) Units - decimal degrees 

Sets the margin beyond the plotting limits in (1) from which data will 
be read (i, e. the outer window). This is to ensure adequate data is available for 
the prediction of grid values on the map borders. 

(iii) NY, NX (INTEGER) 

Denotes the number of grid intervals in the North-South and East-West 
directions (respectively). Values which produce a 1° grid spacing appear to be 
optimum for areas well covered' by altimetry data. 

Note : The arrays X, Y, Z and SD must be dimensioned NY+ 1, 

NX + 1 in the main program and in subroutines AXPLT and PRDT. 

(iv) lAO (INTEGER) 

If lAO = n, every n th point within the inner windown will be plotted 
and identified. 

If lAO = 0, neither points nor tracks will be plotted. 

If lAO 7000, the tracks will be plotted in dashed lines, but no points 
will be plotted. 

(V) AXLAB (REAL*A) Deciznal degrees 

Spacing for labelling of axes. For example, if AXLAB = 5.0, a tick 
and a latitude or longitude value will be printed every 5 degrees. 
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(Vi) IND (INTEGER) 


This enables selection of the projection for the plot (see Program MAIN 

for details). 

(vii) INDl (INTEGER) 

Shoreline plot indicator. If INDl = 0, shoreline will not be plotted. 

If INDl >0, shoreline w.ill be plotted. 

Card 2. Unformatted 

(i) SC FAC (REAL*4) 

Scale f^tor applied to the plot to reduce it to appropriate dimensions. 
This will depend on the extent of the plot (note (i), Card 1), &e chosen projection 
and the desired size of the plot. If the scale chosen is, for example, 1/250000 
then SC FAC = 250Q00.0. 

(ii) CNTINT (REAL*4) 

Contour interval in meters. This should be at least twice the expected 
accuracy of the data. 

(iii) NLC, NDL, HCL, SL, KSM and WTCP. Refer to subroutine CGTOUR. 


A. 2. Output 


(i) The program prints out the details read on cards 1 and 2. It also lists in 
inches the limits of the' plot, and the increments of the grid in N-S and E-W 
directions. 

(ii) The data accepted for use in the prediction (i. e. data accepted into the outer 
window) is now listed. The sequence of the listing is; (i) sequence number within 
window, (ii) observation identifying number, (iii) latitude and (iv) longitude in 
decimal degrees, (v) undulation value, (vi) easting £ind (vii) northing of point at 
plotting scale in decimal inches and (viii) a flag, indicating if the point is the first 
of a section within an orbit. 

(iii) The initial and final points used to plot segments of the tracks is now listed. 

(iv) Details of data selected for the prediction of each grid point is listed row by 
row, row 1 being the row of least latitude. 
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{V) The predicted values of the undulation and the error of prediction are now 
summarized in table form for convenience. 

(vi) Finally, the details of data used for plotting the a?KS are now listed. 


A. 3. Check List When Using the Contouring Program 


Check 

(i) 2nd last data card to certify (a) Map Limits, (b) No. of Grid Intervals in Y and X, 
(c) Track Plot Indicator lAO and (d) Projection choice are correct, 

(ii) DIMENSIONS of X, Y, Z and SD in MAIN program and subroutines PRDT and 
AXPLT, The dimensions of these arrays should be NY+ 1, NX+ 1, 

(iii) DATA card in MAIN program; INC should be set so that no. of data allowed 
into outer window does not exceed 7,000. Note if INC = n, only every n th point 
in window is used, 

(iv) Tape and Tape File Numbers are correct for the area to be mapped. 
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Appendix B; Flow Chart -Main Program 
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